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ABSTRACT 

An  MNDO  study  of  boron-nitrogen  analogues  of  buckminsterfullerene  is  presented.  The 
relative  properties  of  ({2C^q)  ,  ((aB2Ccg)  ,  ((aN2C5g;  .  ^^BNCjg)  ,  23241^2^)  and 

are  studied.  The  heats  of  formation  of  such  60-acom  systems  from  benzene, 
naphthalene  and  their  BN  analogues  are  compared.  It  is  found  that  all  these  hybrids 
ire  approximately  as  stable  as  buckminsterfullerene.  Surprisingly,  it  is  predicted 
-hat  ((SBgQiNgQ)  will  be  stable  and  should  be  relatively  simple  to  synthesize  from 
horazine . 
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I .  INTRODUCTION 

Recent  experiments  iiave  snown  tnat  individual  boron  atoms  tan  displace  carbon 

atoms  in  the  buckv  bail  molecule,  creating  "dopey  ball"  fuilerene  structures  with 

the  stoichiometry  ^  ranges  from  1  to  6,^  The  '3  symbol  indicates  a 

closed  fuilerene  structure  with  all  atoms  forming  an  integral  part  of  the  cage 
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framework.  While  it  has  not  yet  been  demonstrated,  it  seems  reasonable  chat  BN- 
analogues  of  buckmlnsterfullerene  could  also  be  synthesized.  It  has  been  suggested^ 
that  ((SB^^qN^q)  would  not  be  stable  since  it  would  require  Chat  ^3-N  and  B-B  bonds 
exist  in  the  molecule,  presumably  destabilizing  it.  In  particular,  it  is  supposed 
that  N-N  o-bonds  would  be  unfavorable. 

In  response  to  this  suggestion,  we  have  proposed  a  molecule  in  which  Che 
largest  possible  number  of  carbons  are  substituted  subject  to  Che  constraint  that 
there  be  no  B-B  or  N-N  bonds . ^  This  turns  out  to  imply  that  each  pentagon  must 
contain  one  carbon  atom,  yielding  the  stoichiometry  (@0^2224^24^'  optimal 

structure  is  shown  in  Fig.  1,  and  belongs  to  the  point  group.  .A.  simple  Hiickei 
alcuiacion  indicated  that  it  would  be  more  stable  chan  buckmlnsterfullerene. 

.'iothing  in  what  follows  changes  chat  conclusion,  although  we  now  have  considerably 
more  information  to  report. 

We  describe  here  the  results  of  an  MNDO  study  of  boron-nitrogen  analogues  of 
buckmlnsterfullerene.  Given  computational  limitations,  we  are  essentially  limited 
to  closed  shell  systems,  and  hence  are  restricted  to  an  even  nunber  of  substituted 
atoms.  Thus  we  consider  ((3C^q)  ,  ((3B2C5g),  (@N2C,^g)  ,  ((3BNC^g),  ({aCj^2®2^^’24^ 

((SB^qN^q) .  There  are  obviously  many  possible  isomers  of  most  oi  these  species,  and 
we  consider  several  of  each. 

In  order  to  establish  the  accuracy  of  our  method,  we  have  compared  the  MNDO 
calculated  heats  of  formation  for  benzene,  naphthalene  and  borazine  with 
experimental  values.  It  is  seen  that  the  MNDO  calculation  overestimates  all  of 
these  formation  enthalpies  by  about  7X.  The  data  are  given  in  Table  I,  along  with 
those  for  two  naphthalene  analogues.  While  we  have  found  no  comparable  experimental 
data  for  the  latter,  this  information  will  be  useful  in  what  follows. 

The  first  MNDO  calculation  for  buckmlnsterfullerene  was  performed  by  Newton 
and  Stanton.^  Their  value  for  the  standard  heat  of  formation  agrees  with  ours,  +869 
kcal/mol.  A  6-31G  SCF  calculation^  yielded  a  value  of  +672  kcal/mol,  whereas 
classical  methods  gave  +574  kcal/mol  (MM3)®  and  +286  kcal/mol  (MMP2).^  Beckhaus  eC 
al  report  the  value  measured  from  bomb  calorimetry  of  crystalline 
buckmlnsterfullerene  as  +545  kcal/mol.^®  Thus  it  appears  that  MNDO  significantly 
overstates  the  heat  of  formation.  Given  the  relative  accuracy  of  the  method  for 
simple  rings,  one  presumes  that  MNDO  overestimates  the  strain  energy  involved  in 
sphere  closure.  If  this  is  true,  then  our  results  should  nevertheless  be 
qualitatively  correct  since  the  strain  energy  among  our  various  isomers  should  be 
relatively  constant.  On  the  other  hand,  classical  force -field  methods  probably  do 
not  adequately  measure  differences  in  conjugation  energy,  and  hence  are  less  useful 
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for  comparing  molecules  of  similar  geometry.  In  aodition,  they  give  no  Information 
about  charge  distribution. 

The  next  section  contains  the  results  and  discussion  of  our  work.  Section  111 
is  a  brief  conclusion. 

II.  RESULTS  AND  DISCUSSION 

We  first  report  the  relative  stabilitv  of  various  species.  These  data  are 
summarized  in  Table  II,  where  calculated  enthalpies  of  the  following  reactions  are 
compiled: 


10  ---- 

->  (@CgQ)  +  30  H3 

(la) 

10  B3N3Hg  -- 

--->  (@B3qN3q)  +  30  H2 

(lb) 

6  C^gHg  ---- 

->  (@CgQ)  +  24  H, 

(2a) 

6  BjNgHg  --- 

<®»3CI''30>  *  “ 

(2b) 

6  C2B4N^Hg  - 

---->  (@C^232,.M2,^)  f  24 

(2c) 

Each  of  these  reactions  corresponds  to  the  synthesis  of  a  60-acom  fullerene  from  its 
benzene-  and/or  naphthalene-like  precursor.  Among  the  reactants,  benzene, 
naphthalene  and  borazine  are  well  known,  lii*.  reactant  in  Eq .  2b  is 

H  w 

111  '■) 

N  K 

'  » 

1  I. 

and  its  synthesis  from  the  pyrolysis  of  borazine  has  been  reported. The  reactant 
in  Eq.  2c  is 


1  <  1 


(ID 


It 


It 


e  similar  has  been 


This  precursor  has  not  been  made,  thougn  somechinz  quit 
12 

reported. 

From  Table  II  we  see  chat  our  original  conjecture'  that  52^1124^  (which  we 

nickname  CBN-ball)  Is  stable  is  supported  by  the  MNDO  calculations.  As  described  in 
Ref.  3  and  shown  in  Fig.  1,  this  structure  has  one  carbon  on  each  pentagon,  with  two 
pentagons  linked  by  a  C-C  bond.  There  is  symmetry,  with  one  C,  axis.  Thus  one 
can  distinguish  between  polar  and  equatorial  carbon  atoms,  there  being  six  of  each. 

Most  important  are  the  data  for  the  synthesis  of  CBN-ball  from  the  naphthalene 
derivatives,  shown  in  Table  III:  +639  kcal/mole  for  bucky  ball  .’s .  +265  kcal/mole 
for  CBN-ball  (compare  Eqs.  2a  and  2b).  This  corresponds  to  a  difference  of  346 
kcal/mol  in  favor  of  the  substituted  derivative.  Part  of  the  reason  for  this  result 
appears  to  be  the  relative  instability  of  that  analogue  (I),  and  suggests  the 
possibility  of  successful  synthesis  via  the  naphthalene  precursor.  For  comparison, 
we  indicate  the  enthalpy  of  the  reaction  of  bucky  ball  from  benzsne,  and  it  is  seen 
in  Table  III  to  be  17  kcal/mole  nigher  than  the  formation  of  C,-^  from  naphthalene. 
Vet  CgQ  can  be  made  from  the  pyrolvsis  of  benzene  and  other  hydr ocarbons . ‘  There 
is  no  benzene  analogue  with  the  same  stoichiometrv  as  CBN-ball. 

Table  III  also  gives  the  enthalpy  of  reaction  for  the  syntriesis  of 
(denoted  BN-ball)  from  borazine  and  the  naphthalene  analogue.  This  value  is  for  the 
formation  of  an  Isomer  (denoted  A,  shown  in  Fig.  2)  which  is  approximately  as  stable 
as  bucky  ball.  It  contains  6  B-B  and  6  N-N  bonds,  and  is  generated  by  substituting 
all  the  C-C  bonds  in  CBN-ball  with  3-N  moieties.  This  can  be  done  in  a  number  of 
different  ways,  for  example  by  placing  all  the  borons  in  polar  positions  and  the 
nitrogens  along  the  equator.  Ue  have  not  checked  all  possibilities,  but  suppose 
that  they  will  not  differ  significantly  in  energy.  They  may  be  quite  different  in 
their  reaction  chemistries,  however.  The  original  hypothesis,^  that  N-N  and  B-B 
bonds  will  tend  to  destabilize  the  structure,  is  supported  by  our  calculation. 
Results  for  a  second,  less  stable  Isomer  (denoted  B,  with  C^j^symmetry)  are  also 
given  in  Tables  II  and  III.  B  is  obtained  by  replacing  carbons  with  nitrogens  in 
the  northern  hemisphere,  and  with  borons  in  the  southern  hemisphere.  B  is  nearly 
200  kcal/mol  less  stable  than  A,  as  it  contains  9  N-N  bonds  vs.  6  N-N  bonds  for  A. 
(The  number  of  B-B  bonds  in  all  cases  equals  the  number  of  N-N  bonds.)  Ue  also 
checked  intermediate  cases  where  some  C-C  bonds  were  replaced  by  B-N  moieties,  and 
others  by  B-B  and  N-N  moieties,  thus  producing  structures  with  8  or  7  N-N  bonds. 
These  results  confirm  the  notion  that  N-N  bonds  are  destabilizing.  Presumably, 
(^B^qNjq)  could  be  made  more  stable  by  substituting  carbon  atoms  for  each  B-B  and  N- 
N  bond,  until  the  optimal  stoichiometry  of  CBN-ball  is  reached. 

We  now  turn  our  attention  to  the  less  substituted  species,  (@C^gB2),  (@CjgN2) 
and  (^CjgBN) .  Both  of  the  nitrogen-containing  compounds  have  lower  total  energies 
than  bucky  ball,  whereas  ((ac^gBj)  has  higher  energy.  (Note  that  heats  of  formation 
are  not  directly  comparable,  since  the  reference  points  are  not  the  same.)  This  is 
easily  explained  by  a  comparison  of  atomic  sizes  and  corresponding 


eieccronegacivi  ties  .  ;'iLrcp.en  is  .note  e  iec  crcneeac  ive  chan  eich^r  carbon  or  boron, 
ind  hence  its  total  electronic  energy  is  lower.  This  is  not  tiec issar i ly  a  good 
indication  of  the  relative  staoility  of  tr.e  cluster.  .::nall./y  <'•  r  ar^"*  have  reported 
that  the  attempted  synthesis  of  N-containing  species  from  N-doped  graphite  produces 
gas  rather  than  substituted  bucKV  ball.  He  suggests,  and  our  calculations 
confirm,  that  beginning  from  carbon-containing  nitrogen  compounds  will  probably  be 
more  successful.  For  example,  pyrolysis  of  benzene  doped  with  p/ridine  might 
produce  the  desired  result.  However,  verv  recent  work  has  .shown  that  the 
-apor ization  of  graphite  under  T.itrogen  or  ammonia  gas  vield  nitrogen  adducts  of 
buckv  ball,  with  evidence  chat  some  are  substituted  into  ttie  cags,*" 

For  computer-time  reasons,  cur  calculation  is  restricted  to  closed-shell 
-ivstems.  Hence  we  performed  calculations  onlv  wnen  an  even  rumbsr  of  carbons  are 
substituted.  Tur  intention  was  to  probe  a  variety  of  different  isomers,  and  since 
the  available  experimental  results  seem  to  Indicate  that  boron  dopant  atoms  are  non- 
idiacent  on  the  buckv  hall  surface.^  ve  wanted  to  e.stimate  the  eoergy  difference  as 
'.  function  of  dooant  atom  separation  oy  scudving  different  isomers  of  (i^C^gb^j. 

"’nis  has  proved  not  to  be  possible  since  widelv  -paced  dooant  atoms  vields  an  open- 
shell  diradlcai.  This  indicates  that,  like  polymers,  liooco  buck/  ball  has  a 
ieiocallzation  length  beyond  which  electronic  defects  .are  not  smoothed  out.^^  (This 
does  not  contradict  the  well-known  result  that  the  electronic  structure  of  buckv 
bail  is  globally  delocalized,  since  the  diradlcai  requires  electron  deficient 
centers.)  '^Jhiile  this  presents  a  serious  problem  for  our  calculation,  it  would  seem 
to  be  a  boon  for  chemists  interested  in  producing  exo-derivatives  from  these 
ipecies.  They  would  appear  to  be  quite  reactive  and  to  have  an  slectronic  structure 
very  different  from  buckv  ball.  This  is  confirmed  by  a  recent  Car-Parrinello  study 
of  (@C3gB)  and  ((aC^gN).^^ 

For  systems  isoelectronic  with  bucky  ball,  e.g.,  ((ac^gBN) ,  we  were  able  to 

place  the  two  dopants  anywhere  in  the  molecule.  Three  cases  were  studied;  (1)  a 

3-N  bond  between  two  hexagons;  (2)  the  B  and  N  atoms  not  bonded  out  on  the  same 

hexagon;  and  (3)  the  B  and  N  atoms  on  opposite  sides  of  the  molecule.  In  the  first 

case,  the  boron  was  found  to  be  electron  deficient  and  the  nitrogen  electron  rich. 

This  is  consistent  with  the  relative  e Lectronegati/lties ,  and  is  also  the  same  as 

the  pattern  universally  observed  for  all  BN-ball  isomers.  There  is  a  general 

tendency  fer  ;r-eiectron  densitv  to  be  transferred  toward  the  nitrogen,  whereas  n- 
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electrons  move  toward  the  boron.  For  boron,  the  charge  density  in  (OC^gBN)  is 
+0.0523,  and  for  nitrogen,  it  is  -0.1450.  This  compares  with  th=  charge 
distributions  in  BN-ball,  which  range  around  ±  0.23,  though  because  of  the  asymmetry 
there  Is  considerable  variation. 

If  the  boron  and  nitrogen  atoms  are  not  bonded  to  one  anotner,  however,  then 
there  is  obviously  little  o-electron  transfer  between  them.  Nitrogen  is  still 
electronegative,  and  so  borrows  electrons  from  surrounding  carbons.  The  boron  forms 
a  hole  into  which  n-electrons  sink,  and  therefore  becomes  negatively  charged  despite 


its  relatively  low  e  lectronegativitv .  For  th.e  B,  on  same  hexacon.  the  charges  are 
-0.1467  and  -0.0491  respectively,  whereas  tor  atoms  on  r)pposice  .sides  of  the  cluster 
the  ccmparable  values  are  -0.L1?2  and  -0  0942.  In  botti  cases  tl.e  boron  is  more 
negative  than  the  nitrogen.  The  energies  of  the  three  isomers  of  ((3Cc;gBN)  are 
almost  identical,  with  the  B-N  bonded  structure  slightly  favored. 

III.  CONCLUSION 

while  we  reiterate  here  the  semi-quantitative  nature  of  our  results,  it  must 
nevertheless  be  stated  that  !!ND0  Is  best  suited  Co  our  present  purpose.  Our 
Intention  is  to  compare  a  large  number  of  different  structures  where  electrons  are 
globally  delocalized.  This  appears  to  preclude  classical  calculations,  which  would 
not  properly  account  for  the  delocalization.  Assuming  chat  the  discrepancy  between 
MNDO  and  experiment  accrues  mostly  from  strain  energy,  and  noting  that  the  strain 
energy  must  be  approximately  constant  among  all  structures  studied,  MNDO  would 
appear  to  be  Che  method  of  choice. 

The  essential  qualitative  result  mav  thus  oe  •-rummar ized  as  follows:  It  seems 
that  BN  -  subsci  tuted  tierivacives  nf  will  be  .stable,  t-speciallv  those  .species 

which  are  isoeleccronic  with  bucky  oall.  From  the  reactions  described  by  Eqs .  2a-c, 
It  seems  that  (0C|^2^24‘^24 ^  more  stable  than  buckv  ball,  which  is  slightly  more 
stable  chan  (^B^qN^q) .  The  significantly  lower  reaction  enthalpy  for  the  formation 
of  CBN-ball  is  due  in  part  to  the  relative  instability  of  the  naphthalene- like 
precursor.  But,  based  on  total  energies  (Table  II).  ail  species  rival  bucky  ball 
with  respect  to  thermodynamic  stability. 

Kinetic  stability  is  harder  to  judge.  As  was  noted  above,  lucky  ball  has  been 

synthesized  by  the  pyrolysis  of  benzene.  ^  This  suggests  that  ((JB^qN-^q)  might  be 

obtained  from  borazine  either  by  thermal  degradation  as  per  Eq .  lb,  or  by  high- 

19 

temperature  air  oxidation.  ^@^12^24^24^  represents  a  more  difficult  synthetic 

problem.  Precursors  having  C;B:N  atomic  ratios  of  1:2:2  would  be  especially 

attractive,  as  Illustrated  by  Eq .  2c,  but  such  species  are  rare.  However,  a 

12 

synthetic  route  to  compound j  having  ftised  CB2N2  rings  has  been  described. 

Once  made,  both  BN-ball  and  CBN-ball  are  expected  to  have  significantly 
different  chemical  properties  from  those  of  (^C^q) .  Bucky  ball  l:self  is  reactive 
toward  both  nucleoDhiles  and  electrophiles.'^^  The  uneven  charge  distributions  of 
the  B-  and  N-doped  clusters  should  Increase  both  types  of  reactivity,  making  them 
attractive  as  ligands,  electron  transfer  agents,  etc.  They  would  be  important 
additions  to  the  chemical  arsenal. 
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Table  I;  Calculated  vs.  Experimental  Enthalpies  of  Formation  for  Selected 


Compounds . 

Compound 

(kcal/mol) 

Reference 

MNDO 

Expt. 

^6^6 

-  21.3 

19.8 

*-38.3 

-*-56.1 

- 

33N3H6 

-131.1 

-124 

5 

-132 

- 

33N3H3 

■  .  i. 

- 

*"60 

*-869.3 

‘545 

10 

Table 


II:  Calculated  Total  Enercies  and  Enthalpies  of  Formation  for  Some  Sixtv-Atom 

Clusters . 


Compound 

Geome  trv 

Energy  ',eV) 

(kcal/ 

-  76 36 

^869.3 

^60 

'h 

®30‘'30 

■S  N-N  bonds 

-8360 

-681 . = 

^  'I 

^30^  30 

3  M - N  bonds 

-3359 

-  5  5 1 .  7 

®30''30 

9  N-N  bonds 

-8351 

-989,0 

r 

"s 

-8221 

-  526  .  8 

'^3  s'"? 

-p  hono 

-  7795 

♦  888,  } 

C.gTo 

-h  bond 

-  7  795 

1-Q05  2 

CggBj 

h-h  bond 

-  7521 

8  6  9  9 

CggBN 

h-h  bond 

-7660 

f839  9 

CggBN 

same  hex 

-7659 

-h861 ,6 

CggBN 

opposite 

-  7658 

-t-881. 5 

formation  from 

the  standard 

states  of  carbon, 

boron  and, 

25’’C. 
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able  III: 


Enthalpies  for  oeleccea  Synthetic  Reactions  Producing  Si:<tv-Acom 
Clusters,  as  Calculated  from  MNDO  Data  in  T.ibles  I  .inn  II. 


From  Benzene  or  .Analogue 

10  CgHg  . >  ((ac^Q;  -t-  .30  Ho 

10  B3N3Hg  . >  +  30  H2 

From  Naphthalene  or  Analogues 

^  *^10^8  . "  ^^^60^  ^  ^2 

5  B^NgHg  - >  ((aSg^Ngg)  +  24  H2 

6  . >  '  3r,  _  3o  .  I'o  /  )  ^  24  H„ 

^Isomer  A  of  ((cB-qN-,q)  . 

^Isomer  B  of  (@B.,q;;,q)  . 


aH°  M'.ca  l/mole  ) 

+  656.1 

f62'5..5'^  r822.0^ 

aH°  < /mo  le  ) 

^633. 3 

rhlV.:^  +864.0^ 
--.165. 2 


FIGURE  CAPTION'S 


Figure  1 


Figure  2 


The  scruccure  oC  (;3C,  .  )  as  viewed  along  :.ho  C,  rotational  axis. 

Nine  of  the  twelve  carbons  (all  six  equatorial  plus  three  polar)  are 
shown.  Taken  from  Ref.  3. 

The  A  isomer  of  ((SB^qN^q)  •  Borons  are  shown  as  hatched  circles,  the 
others  being  nitrogen  acorns.  Note  that  if  a  six-memberea  ring  has  an  N-N 
bond,  it  must  also  have  a  B-B  bond  in  order  to  be  i soe lec tronic  with  Che 
carbon  equivalent.  .About  half  the  molecule  is  shown. 


Carbon 


